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Abstract 

Numerical  simulations  of  a  50  cm2  polymer  electrolyte  fuel  cell  (PEFC)  with  36  channels  are  carried  out  to  study  the  lateral  transport 
of  moisture  and  reactant  between  two  neighboring  channels  with  counter  flow  on  the  cathode  side.  Massive  computations  with  2.7  million 
computational  elements  are  performed  to  capture  the  intricate  electrochemical  and  transport  phenomena  in  PEFC  in  three-dimensions.  Two 
cases  are  examined  and  compared.  One  is  under  the  common  assumption  that  molecular  diffusion  dominates  the  species  transport  in  the 
gas  diffusion  layer  (GDL)  and  thus,  the  lateral  convection  effect  is  neglected.  The  other  case  is  to  include  the  flow  and  convection  effects  in 
the  GDL  using  a  realistic  permeability  of  the  porous  GDL.  Numerical  results  elucidate  the  mechanism  and  extent  of  internal  humidification 
induced  by  lateral  moisture  diffusion.  In  addition,  it  is  found  that  given  the  typical  GDL,  there  exists  a  large  pressure  drop  between  two 
adjacent  channels  in  counter  flow,  one  flowing  in  from  the  inlet  and  the  other  flowing  out  to  the  outlet,  causing  severe  reactant  bypass  between 
the  two  flow  paths  of  reverse  direction.  The  bypass  results  in  reactant  flow  “short-circuit”  and  greatly  diminishes  the  internal  humidification 
benefit  intended  by  the  counter-flow  design. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Due  to  development  of  perfluorosulfonic-acid  membranes 
and  the  dramatic  reduction  in  precious  metal  catalyst  loading, 
the  polymer  electrolyte  fuel  cell  (PEFC)  has  become  a  viable 
technology  [1].  A  PEFC  consists  of  gas  channels,  gas  diffu¬ 
sion  layers  (GDLs),  and  a  proton-conductive  membrane  with 
platinum  catalyst  coated  on  the  two  surfaces,  called  the  mem¬ 
brane  electrode  assembly  (ME A).  Gas  channels  are  grooved 
in  graphite  or  metal  plates,  where  injected  reactant  streams 
are  distributed  for  electrochemical  reactions.  Protons  are  pro¬ 
duced  from  hydrogen  oxidation  in  the  anode  catalyst  layer, 
and  pass  through  the  membrane,  carrying  water  molecules 
via  electro-osmotic  drag,  to  the  cathode  catalyst  layer  where 
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the  oxygen  reduction  reaction  (ORR)  occurs,  with  water  as 
the  product. 

Water  management  is  a  central  issue  in  PEFC  technol¬ 
ogy  because  high  water  content  results  in  high  proton  con¬ 
ductivity  of  the  membrane,  thus  decreasing  ohmic  voltage 
loss.  Commonly,  two  methods  are  used  to  hydrate  the  mem¬ 
brane.  One  is  to  externally  humidify  inlet  reactant  gases.  This 
strategy  requires  a  humidifier,  making  it  complex  and  ex¬ 
pensive  to  be  applied  to  portable  or  mobile  systems.  The 
other  way  is  to  make  full  use  of  water  production  inside  a 
PEFC  and  humidify  reactant  gases  internally.  In  the  latter  ap¬ 
proach,  Qi  and  Kaufman  [2]  proposed  a  counter-flow  cathode 
flowfield  design  where  an  inlet  gas  channel  is  placed  next 
to  an  outlet  channel  in  hopes  that  the  moist  air  in  the  out¬ 
let  channel  will  help  humidify  the  dry  air  in  the  inlet  chan¬ 
nel.  Such  a  self-humidification  process  is  of  technological 
interest  as  it  has  the  potential  to  feed  dry  reactants  without 
external  humidification.  While  cell  polarization  experiments 
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Nomenclature 

a 

effective  catalyst  area  per  unit  volume  (m2/m3) 

A 

superficial  electrode  area  (m2) 

Ck 

molar  concentration  of  species  k  (mol/m3) 

D 

mass  diffusivity  of  species  (m2/s) 

F 

Faraday’s  constant,  96,487  C/equivalent 

io 

exchange  current  density  (A/m2) 

I 

current  density  (A/cm2) 

j 

transfer  current  (A/cm3) 

K 

permeability  (m2) 

n 

the  direction  normal  to  the  surface 

nd 

electro-osmotic  drag  coefficient  (H20/H+) 

P 

pressure  (Pa) 

R 

gas  constant  (8.134  J/mol  K) 

s 

stoichiometric  coefficient  in  electrochemical 
reaction 

S 

source  term  in  transport  equations 

t 

time  (s) 

T 

temperature  (K) 

ll 

velocity  vector  (m/s) 

u0 

open  circuit  potential  (V) 

Well 

cell  potential  (V) 

mole  fraction 

Greek 

a 

transfer  coefficient  or  net  water  transport  co¬ 
efficient  per  proton 

£ 

porosity 

0 

phase  potential  (V) 

T) 

surface  overpotential  (V) 

K 

ionic  conductivity  (S/m) 

X 

membrane  water  content 

[l 

viscosity  (kg/m  s) 

a 

electronic  conductivity  (S/cm) 

P 

density  (kg/m3) 

X 

shear  stress  (N/m2,  or  tortuosity  factor) 

h 

stoichiometric  flow  ratio 

Superscripts  and  Subscripts 

a 

anode 

avg 

average 

c 

cathode 

e 

electrolyte 

eff 

effective  value 

g 

gas  phase 

in 

inlet 

k 

species 

m 

membrane 

0 

standard  condition,  273. 15  K  and  101.3  kPa 
(1  atm) 

ref 

reference 

s 

solid  phase  of  electrode 

sat 

saturate  value 

w 

water 

have  been  reported  by  Qi  and  Kaufman  [2]  for  this  flowfield 
design,  a  fundamental  understanding  of  the  complex  trans¬ 
port  phenomena  uniquely  present  in  this  flow  configuration  is 
absent. 

Mathematic  models,  describing  detailed  physico¬ 
chemical  processes  and  coupling  transport  equations  with 
electrochemical  kinetics,  are  an  important  tool  to  aid  in 
fundamental  understanding  of  relevant  phenomena  in 
PEFCs.  To  date,  mathematic  models  have  been  developed 
from  the  simple  one-dimensional  [3-5],  two-dimensional 
[6,7]  to  the  complex  three-dimensional  models.  Springer  et 
al.  [4,5]  presented  a  detailed  description  of  water  behavior 
in  the  membrane  and  experimentally  measured  the  water- 
diffusion  coefficient,  electro-osmotic  drag  coefficient  and 
membrane  proton  conductivity  as  functions  of  membrane 
water  content.  Fitted  curves  were  incorporated  into  their  one¬ 
dimensional  model  to  capture  water  transport  phenomena  in 
the  membrane  via  molecular  diffusion  and  electro- osmotic 
drag.  Fuller  and  Newman  [6]  considered  two-dimensional 
MEA  in  the  through-plane  and  flow  directions  operat¬ 
ing  on  reformed  methanol.  Water  management,  thermal 
management  and  utilization  of  fuel  were  examined  in 
detail. 

A  comprehensive  non-isothermal,  three-dimensional 
model  was  developed  by  Berning  et  al.  [8].  They  simulated  a 
half  of  a  single  channel  to  study  the  temperature  distribution. 
However,  their  model  ignored  the  membrane-electrode  as¬ 
sembly,  thus  not  capturing  the  transport  and  electrochemical 
mechanisms  in  these  vital  components  of  PEFCs.  The  same 
assumption  also  existed  in  the  model  of  Dutta  et  al.  [9,10]. 
An  important  contribution  of  Dutta’ s  work  is  to  simulate  a 
fuel  cell  with  20-branch  serpentine  channels.  However,  only 
about  200,000  grid  points  were  employed  in  the  simulations, 
raising  the  concern  of  grid  dependence  and  numerical  inac¬ 
curacy.  In  addition,  the  governing  equations  were  employed 
only  in  the  gas  channels  and  GDFs,  with  the  MEA  excluded. 
Water  transport  in  the  MEA  was  analytically  modeled  based 
on  the  approximations  of  constant  transport  properties,  and 
the  diffusion  flux  through  the  membrane  was  calculated  based 
on  the  difference  of  water  concentrations  in  the  two  catalyst 
layers.  In  contrast,  Wang  and  co-workers  [11-15]  presented  a 
single-domain  model,  including  detailed  mechanisms  in  the 
MEA,  such  as  electrolyte  potential,  current  and  water  distri¬ 
butions  in  the  membrane  as  well  as  catalyst  layers.  Thus,  their 
model  is  chosen  for  this  study  of  internal  humidification  and 
membrane  hydration. 

The  goal  of  this  paper  is  to  investigate  the  internal 
humidification  phenomenon  through  the  simulations  of  a 
50  cm2  PEFC  with  counter  flow  in  the  cathode  flowfield. 
The  geometry  of  the  fuel  cell  is  similar  to  the  one  in 
the  work  by  Qi  and  Kaufman  [2]  except  that  a  complete 
PEFC  with  36  channels  is  considered.  We  expand  upon 
the  model  of  Um  et  al.  [11]  and  numerically  implement 
it  into  Star-CD®,  a  commercial  computational  fluid  dy¬ 
namics  (CFD)  software  package,  based  on  its  user  coding 
capability. 
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2.  Numerical  Model 

Fig.  1(a)  shows  the  computational  domain  of  the  double¬ 
path  PEFC  with  36  channels  under  consideration.  There  are 
two  parallel  branches  of  18  channels  in  a  serpentine  fash¬ 
ion  on  each  side.  Anode-side  and  cathode-side  flow  channels 
are  symmetric  and  properly  aligned  on  top  of  each  other. 
Fig.  1(b)  shows  part  of  the  mesh  near  inlet  and  outlet  ar¬ 
eas.  Finer  grids  are  used  to  capture  the  transport  phenomena 
in  channel  U-turns  and  inlet  regions.  Detailed  mesh  in  the 
through-plane  direction  of  the  MEA  is  plotted  in  Fig.  1(c). 
Three-dimensional  mesh  is  necessary  to  capture  the  impor¬ 
tant  electrochemical  and  transport  mechanisms  occurring  in 
the  MEA  and  GDL.  These  include  species  transport  in  the 
through-plane  direction,  flow  and  reactant  consumption  in 
the  flow  direction,  and  gas  flow  and  water  transport  in  the 
in-plane  direction  in  the  GDL  over  the  land  separating  two 
neighboring  channels.  Fig.  1(d)  schematically  shows  the  top 
view  of  the  computational  domain  of  double-path  PEFC  with 
36  channels  and  the  positions  of  inlets  and  outlets.  In  order  to 
focus  on  the  effects  of  internal  humidification,  we  consider 
counter  flow  only  in  the  cathode  flowfield  and  use  a  common 
co-flow  pattern  in  the  anode  one.  In  addition,  dry  air  is  fed 
into  Channels  1  and  35  on  the  cathode,  and  fully-humidified 
hydrogen  is  injected  into  Channels  1  and  2  for  the  anode.  The 
mid-length  cross-section  of  the  fuel  cell  is  typically  selected 
to  show  the  detailed  profiles  of  flow  and  transport  parame¬ 
ters.  To  distinguish  the  two  flow  paths  in  the  cathode  side, 
the  one  with  Channel  1  as  the  inlet  is  defined  as  Path  1  while 
the  other  one  starting  from  the  other  end  is  Path  2.  The  PEFC 
geometry  is  shown  in  Table  1 . 


Table  1 


Geometrical  and  physical  parameters 


Quantity 

Value 

Gas  channel  depth  (mm) 

1.0 

Gas  channel  width  (mm) 

1.0 

Shoulder  width  (mm) 

1.0 

Diffusion  layer  thickness  (mm) 

0.3 

Catalyst  layer  thickness  (mm) 

0.01 

Membrane  (N112)  thickness  (mm) 

0.051 

Fuel  cell  height/length  (mm) 

72.0/72.0 

Anode/cathode  pressure,  p  (atm) 

2. 0/2.0 

Reference  current  density,  7ref  (A/cm2) 

1.0 

Stoichiometric  ratio  £  in  the  anode/cathode 

2. 0/2.0 

RH  of  anode/cathode  inlet  (%) 

100/0 

Temperature  of  fuel  cell,  T  (K) 

353 

Porosity  of  the  GDL,  s  [18] 

0.6 

Porosity  of  catalyst  layers,  eg  [18] 

0.4 

Volume  fraction  of  ionomer  in  catalyst  layers, 

0.26 

Gn  [18] 

Permeability  of  the  GDL,  K  (m2)[18] 

10“12 

H2  diffusivity  in  anode  gas  at  standard  condition, 

1.1028  x  10“4 

A),H2,a  (m2/s)  [22] 

H2O  diffusivity  in  anode  gas  at  standard  condition, 

1.1028  x  10“4 

£>o,w,a  (m2/s)  [22] 

O2  diffusivity  in  cathode  gas  at  standard  condition, 

3.2348  x  10“5 

D0,o2,c  (m2/s)  [22] 

H2O  diffusivity  in  cathode  gas  at  standard  condition, 

7.35  x  10“5 

£>o,w,c  (m2/s)  [22] 

2.7.  Model  equations 

The  single-phase  PEFC  model  consists  of  non-linear,  cou¬ 
pled  partial  differential  equations  describing  the  conservation 
of  mass,  momentum,  species  and  charge  with  electrochem- 
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Fig.  1.  (a)-(d)  Geometry,  grids  and  top  view  of  the  computational  domain  of  the  double-path  PEFC  with  36  channels. 
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ical  reactions.  The  equations  can  be  written,  in  the  vector 
form,  as  [11]: 

Continuity  conservation  :  V  •  ( pu )  =  0  (1) 

Momentum  conservation  : 

1 

-r- V  •  ( puu )  =  -V p  +  V  •  r  +  Su  (2) 

£z 

Species  conservation  :  V  •  ( uC k)  =  V  •  (DjffVCk)  +  S k 

(3) 

Charge  conservation  :  V  •  (/ceffV0e)  +  Se  =  0  (4) 

where  p ,  u,  p,  Cp  and  0e,  respectively,  denote  the  density, 
superficial  fluid  velocity  vector,  pressure,  mole  fraction  of 
species  k  and  electrolyte  potential.  The  species  considered  are 
hydrogen,  oxygen  and  water.  In  Eqs.  (1)  and  (2),  use  has  been 
made  of  a  constant-flow  assumption  in  which  the  fluid  density 
is  constant  and  the  mass  source/sink  in  the  continuity  equation 


Table  2 


Source  terms  for  the  conservation  equations  in  each  region 


Su 

s k 

Sep 

Gas  channels 

0 

p  -+ 

- u 

K 

0 

— 

Diffusion  layers 

Su  = 

0 

0 

Catalyst  layers 

Su  = 

P  -+ 

- u 

K 

II 

< 

(7'-) 

Skj 

nF 

S<t>=j 

Membrane 

Su  = 

P  - 

- u 

K 

0 

0 

Note:  nd  is  the  electro-osmotic  drag  coefficient  for  water.  For  H2  and  O2, 
nc\  =  0. 


is  ignored.  This  constant-flow  approximation  has  proved  to 
be  valid  under  common  PEFC  operations  [16].  Source  terms, 
electrochemical  properties  and  physical  properties,  identified 
for  various  regions  of  a  fuel  cell,  are  listed  in  Tables  2-4, 
respectively. 

In  the  GDL,  a  small  value  of  permeability  K  will  result 
in  a  large  source  term  Su  in  Eq.  (2),  turning  Eq.  (2)  into  the 


Table  3 

Electrochemical  properties 


Description 


Anode 


Cathode 


Transfer  current  density,  j  (A/m3) 

Surface  overpotential,  ij  (V) 

Equilibrium  potential,  U0  (V) 

Exchange  current  density  x  reaction  surface  area,  aio  (A/m3) 
Transfer  coefficient,  a 


6^0, a 


Cho 


1/2 


Ot  a  +  ac 


Cfb.ref  /  V  RT 

0s  -  0e  -  U0  (with  0S  =  0) 

0 

1.0  X  109 
aa  +  ac  =  2 


F  •  rj 


ai  o,c 


_(  “rt-'i) 

Co2,ref 

0s  “  0c  -  U0  (with  0s  =  Vcell) 

1.23- 0.9  X  10“3  (T- 298) 

10000 

ac  =  1 


Table  4 

Physical  and  transport  properties 

Quantity 

Value 

Water  activity,  a  [4] 


Water  content  in  membrane,  k  [4] 

Ionic  conductivity  of  membrane,  k  (S/m)  [4] 


H2O  diffusivity  in  membrane,  D (m2/s)  [19] 


Electro-osmotic  drag  coefficient,  nd  [20] 


Membrane  density,  p  (kg/m3)  [21] 


CWRT 


a  = 


,sat 


p 

log  10 pSM  =  -2.1794  +  0.02953  {T  —  273.15)  -  9.1837  x  10"5(r-  273.15)2 
+  1.4454  x  10"7(r-  273.15)3 

0.043  +  17.81a  —  39.85a2  +  36.0a3  for  0  <  a  <  1 

for  1  <  a  <  3 

1  1 


A  '  14  +  \A(a  —  1) 
(0.005 139A-  0.00326)  exp 


1268 


303 


Dm  = 


3.1  x  10“3A(eo  m  -  1)  •  e[-2346/r]  for  0  <  A  <  3 
4.17  x  10_4A.(1  +  161e-A)  •  e[_2346/r^  otherwise 


1.0 


for  k  <  14 


nd  =  {  1.5 
~8~ 


(A  —  14)  +1.0  otherwise 


1.98  +  0.0324*  , 

p  =  -  x  103 

1  +  0.0648* 


H2/O2  diffusivity  in  membrane,  (m2/s)  [3] 

Gas  diffusivity  in  porous  media,  Deff  (m2/s) 

Diffusivity  in  gas,  D  (m2/s)  [22] 

Viscosity  of  anode/cathode  gas,  /x  (kg/m  s)  [23] 


2.59  x  10_6/1.22  x  10“6 
Deff  =  s15D 


[i  =  9.88  x  10 6Xh2  +  1.12  x  10 5AH2o  +  2.01  x  10-5  AN2  +  2.3  x  10 5AO2 
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Darcy’s  law  for  porous  media: 

Vp  =  fu  (5) 

In  addition,  notice  that  the  present  model  is  a  single-phase 
approach  in  which  the  total  water  amount  is  tracked  without 
distinguishing  liquid  water  from  water  vapor.  Such  a  single¬ 
phase  model  is  valid  generally  under  the  condition  that  the 
liquid  volume  fraction  is  low  or  liquid  droplets  are  small  to 
form  a  mist  flow.  This  approach  is  particularly  well  suited 
for  the  present  study  since  only  the  low-humidity  operation 
is  concerned. 


2.2.  Boundary  conditions 

Eqs.  (1)— (4)  form  a  complete  set  of  governing  equations 
with  eight  unknowns:  u ,  p ,  Ch2,  Go2,Ch2o  and  0e-  Their 
corresponding  boundary  conditions  are  specified  as  follows. 


2.2.7.  Flow  inlet  boundaries 

The  inlet  velocity  um  in  a  gas  channel  is  expressed  by  the 
respective  stoichiometric  flow  ratio,  i.e.  §a  or  §c,  defined  at  a 
reference  current  density,  7ref,  as: 


CH2  Pawin,a 


Tref-^ 

IF 


and 


Co2  Pc^  in,c 


/reijd 

4  F 


where  Aa  and  Ac  are  the  flow  cross-sectional  areas  of  the 
anode  and  cathode  gas  channels,  respectively.  The  inlet  mo¬ 
lar  concentrations  are  determined  by  the  inlet  pressure  and 
humidity  according  to  the  ideal  gas  law. 


2.2.2.  Outlet  boundaries 

Fully  developed  or  no-flux  conditions  are  applied: 


2.2.3.  Wall 

No-slip  and  impermeable  velocity  condition  and  no-flux 
condition  are  applied: 


2.3.  Numerical  procedure 


The  conservation  equations  are  solved  by  Star-CD  soft¬ 
ware  with  SIMPLE  algorithm  [17],  using  a  parallel  com¬ 
putational  methodology  for  a  Linux  PC  cluster.  The  source 
terms  and  physical  properties  are  incorporated  in  the  user 
codes.  Cell  voltage  is  prescribed  and  local  current  density 
distribution  is  computed  in  all  simulations  to  be  presented 
in  this  work.  Alternatively,  an  average  current  density  can 
be  prescribed  and  the  cell  voltage  is  computed  if  the  elec¬ 
tron  transport  equation  is  included,  as  shown  in  [15].  After 
grid-independence  study,  about  2.7  million  computational  el¬ 
ements  (44  x  100  x  576)  were  determined  to  be  necessary  to 


capture  fine  details  of  flow  and  transport  phenomena  in  this 
PEFC  configuration.  A  simulation  typically  requires  about 
100  iterations  and  takes  nearly  5h  on  9  nodes  of  1.4  GHz 
AMD  Athelon  Thunderbird  CPU  and  5 12MB  DDR  SDRAM. 


3.  Results  and  discussion 

Two  cases  are  compared  in  our  study,  in  Case  1 ,  we  assume 
there  is  no  velocity  in  the  GDL,  thus  species  transport  is 
dominated  by  diffusion.  The  other  one,  Case  2,  is  to  apply  the 
Darcy’s  law  in  the  GDL  based  on  its  permeability  and  explore 
the  convection  effect  on  internal  humidification  additionally. 

Fig.  2  shows  the  local  velocity  distribution  at  the  mid¬ 
depth  of  the  cathode  gas  channels  for  Case  1 .  The  grid  is  much 
finer  at  the  U-turns,  inlet  and  outlet,  in  order  to  capture  the 
complex  flow  in  these  locations,  such  as  flow  recirculation. 
From  the  inlet  of  Channel  1,  the  gas  velocity  profile  under¬ 
goes  a  change  from  the  uniform  flow  to  fully-developed  flow 
which  can  be  seen  in  Channel  4  before  the  U-turn.  A  similar 
parabolic  profile  velocity  is  displayed  at  Channel  2,  which  is 
the  outlet  channel  of  Path  2.  In  the  gas  channel,  the  convec¬ 
tion  dominates  the  mass  transport,  and  species  is  distributed 
down  the  flow  paths. 

Fig.  3  shows  the  distributions  of  average  O2  and  H2O 
concentrations  in  the  cathode  gas  channel  of  each  flow  path 
for  Case  1.  O2  concentration  is  seen  to  decrease  along  the 
channels  from  the  two  inlets,  due  to  the  oxygen  reduction 
reaction.  The  concentration  difference  between  the  two  flow 
paths  is  maximum  in  the  inlet  and  outlet  areas  and  diminishes 
to  zero  in  the  middle  of  the  cell.  In  addition,  O2  displays  a 
nearly  symmetric  profile  between  the  two  flow  paths,  indicat¬ 
ing  that  nearly  equal  amount  of  O2  is  consumed  along  each 
path. 

In  contrast,  H2O  concentration  increases  along  the  chan¬ 
nel  from  the  inlet  to  outlet  due  to  water  production  of  ORR. 
Furthermore,  H2O  concentration  exhibits  asymmetric  distri- 
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Fig.  2.  Velocity  distribution  at  the  mid-depth  of  the  cathode  gas  channel 
without  convection  in  the  GDL. 
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Fig.  3.  Average  O2/H2O  concentrations  in  the  cathode  gas  channel  for  each 

butions  between  the  two  flow  paths.  This  is  due  to  the  effect  of 
anode  co-flow  pattern  where  humidified  H2  is  fed  into  Chan¬ 
nels  1  and  2  from  same  end  of  the  cell.  In  addition,  there  exists 
great  difference  of  water  concentration  between  the  two  flow 
paths  in  the  inlet  or  outlet  areas,  where  internal  humidifica- 


flow  path  without  convection  in  the  GDL  at  Vceii  =  0.65  V  and  /avg  =  0.91  A/cm2. 

tion  is  expected  to  occur  due  to  lateral  diffusion  through  the 
porous  GDL. 

Details  of  the  lateral  transport  of  H2O  and  O2  between 
two  adjacent  channels  in  the  different  flow  paths  are  shown 
in  Fig.  4(a)  and  (b),  respectively.  The  cross-section  is  taken 


Fig.  4.  (a)  H2O  and  (b)  O2  distributions,  respectively,  in  the  cathode  at  the  mid-length  cross-section  between  Channels  1  and  2  without  convection  in  the  GDL 
at  Vceii  =  0.65  V  and  /av g  =  0.91  A/cm2. 
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Fig.  5.  (a)  H2O  and  (b)  O2  distributions,  respectively,  in  the  cathode  at  the  mid-length  cross-section  between  Channels  5  and  6  without  convection  in  the  GDL 
at  Vceii  =  0.65  V  and  /av g  =  0.91  A/cm2. 


on  the  cathode  side  at  the  mid-length  of  the  fuel  cell  between 
Channels  1  and  2  under  Vcen  =  0.65  V  and  /avg  =  0.91  A/cm2. 
A  clear  gradient  in  the  water  concentration  towards  the  dry 
inlet  channel  (Channel  1)  is  evident.  Fig.  4(b)  shows  O2  con¬ 
centration  is  lowest  in  the  GDL  under  the  land.  This  is  due  to 
the  ORR  under  the  land  and  the  fact  that  the  thin  GDL  limits 
the  O2  supply. 

Fig.  5(a)  and  (b)  display  the  corresponding  H2O  and  O2 
distributions,  respectively,  between  Channels  5  and  6  under 
the  same  conditions  as  in  Fig.  4.  Fig.  5(a)  shows  that  there 
also  exist  dramatic  concentration  differences  between  these 
two  channels  both  for  water  and  oxygen.  Similar  to  Fig.  4(a), 
the  highest  concentration  of  H2O  appears  under  the  land. 
Fig.  5(b)  also  indicates  virtually  no  O2  exchange  between 
the  two  channels. 

The  reason  why  there  is  no  direct  water  exchange  between 
the  two  channels  in  Figs.  4  and  5  is  due  to  local  water  produc¬ 
tion  under  the  land.  If  there  were  no  such  water  source,  water 
vapor  would  diffuse  from  the  humid  channel  to  the  dry  one 
under  the  concentration  difference.  However,  water  source 
under  the  land  raises  the  local  water  concentration  even  be¬ 
yond  the  wet  channel  nearby,  thus  preventing  water  diffusion 
between  the  two  channels.  Nonetheless,  a  larger  portion  of 
product  water  is  laterally  removed  into  the  drier  gas  channel, 
so  there  is  a  limited  extent  of  internal  humidification. 


Fig.  6(a)  and  (b)  show  H2O  and  O2  distributions,  respec¬ 
tively,  in  the  cathode  side  between  Channels  35  and  36.  In 
this  location,  different  to  Figs.  4  and  5,  water  concentration 
under  the  land  is  lower  than  the  humid  Channel  36,  so  a  lat¬ 
eral  water  flux  induced  by  diffusion  occurs  through  the  GDL 
between  the  two  channels.  Simultaneously,  water  produced 
under  the  land  also  enters  the  dry  Channel  35.  This  more  ben¬ 
eficial  effect  of  internal  humidification  occurs  in  the  very  last 
channel  towards  the  exit.  Fig.  6(b)  shows  that  there  is  still  no 
O2  exchange  between  the  two  channels.  Rather,  O2  needed 
for  the  electrochemical  reaction  under  the  land  is  supplied 
from  both  channels. 

The  effect  of  internal  humidification  is  shown  quantita¬ 
tively  in  Fig.  7,  which  indicates  the  water  flux  across  the 
two  surfaces  of  the  cathode  GDL  along  the  mid-length  cross- 
section  for  Case  1.  The  areas  under  the  curves,  i.e.  the  in¬ 
tegrals  of  the  curves,  represent  the  amount  of  water  across 
each  face  of  the  GDL  per  unit  length.  Due  to  no  water  pro¬ 
duction  in  the  GDL,  the  area  between  the  two  curves  across 
the  two  surfaces  of  the  GDL  represents  the  net  water  flux 
along  the  lateral  direction.  It  can  be  seen  that  the  amount  of 
water  into  the  gas  channel  is  almost  twice  that  across  the  in¬ 
terface  between  the  catalyst  layer  and  GDL  under  Channel  1 , 
while  the  water  fluxes  across  the  two  surfaces  of  the  GDL  are 
almost  equal  under  Channel  2.  The  different  amount  of  prod- 
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Fig.  6.  (a)  H2O  and  (b)  O2  distributions,  respectively,  in  the  cathode  at  the  mid-length  cross-section  between  Channels  35  and  36  without  convection  in  the 
GDL  at  Veen  =  0.65  V  and  7avg  =  0.91  A/cm2. 


uct  water  captured  by  the  inlet  Channel  1  and  outlet  Channel 
2  indicates  that  internal  humidification  exists  and  has  a  sig¬ 
nificant  effect  on  humidifying  the  dry  inlet  flow.  The  same 
trend  also  exists  in  the  channels  close  to  the  air  inlet  such 
as  Channels  4,  5,  34  and  35,  and  the  channels  close  to  the 


exhaust  outlet  such  as  Channels  3,  6,  33  and  36.  In  Chan¬ 
nel  36,  there  exists  a  negative  flux  across  the  front  face  of 
the  GDL,  which  means  that  water  diffuses  from  Channel  36 
into  the  GDL  and  is  transported  laterally  through  the  GDL  to 
Channel  35.  Similar  phenomenon  also  appears  in  the  outlet 


Fig.  7.  Water  flux  across  surfaces  of  the  cathode  GDL  at  the  mid-length  cross-section  without  convection  in  the  GDL  under  Vcen  =  0.65  V  and  7av  g  =  0.91  A/cm2. 
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Fig.  8.  Net  water  transport  coefficient  per  proton  at  the  mid-length  cross-section  without  convection  in  the  GDL  at  VCeii  =  0-65  V  and  7av g  =  0.91  A/cm2. 


channels,  such  as  Channels  32  and  33.  However,  the  amount 
of  water  transported  between  the  two  neighboring  channels 
remains  very  small.  In  addition,  the  highly  oscillatory  nature 
of  these  surface  water  fluxes  is  indicative  of  the  complex  wa¬ 
ter  transport  phenomena  occurring  in  the  flow  field  and  thus 
calls  for  a  fine  mesh  to  resolve  sharp  gradients  in  the  lateral 
direction. 

Fig.  8  shows  the  distribution  of  net  water  transport  coef¬ 
ficient  per  proton  at  the  same  cross-section.  The  a-value  re¬ 
flects  the  combined  effect  of  electro-osmotic  drag  and  back 
diffusion  of  water,  with  positive  values  being  water  trans¬ 
ported  from  the  anode  to  cathode.  It  is  highly  dependent  on 
the  location  and  varies  from  —0.5  to  1.5.  The  peak  values 
appear  near  the  cathode  dry  inlet  channels,  Channels  1,  4,  34 


and  35,  while  near  the  cathode  outlet  channels,  such  as  Chan¬ 
nels  2,  3  and  36,  a  is  less  than  unity  because  of  strong  back 
diffusion.  In  addition,  under  the  land,  such  as  the  one  between 
Channels  33  and  34,  a  dramatically  changes  from  a  positive 
to  a  negative  value.  This  can  be  explained  by  water  accumu¬ 
lation  under  lands.  Furthermore,  Fig.  8  clearly  indicates  that 
a  is  highly  dependent  upon  the  local  conditions. 

In  Case  1 ,  the  velocity  in  the  GDL  has  been  ignored  and  the 
species  transport  is  only  by  diffusion.  However,  velocity  will 
be  induced  by  the  pressure  gradient  following  the  Darcy’s  law 
in  porous  media.  Fig.  9  shows  the  gas  pressure  distribution  in 
the  cathode  GDL  at  the  mid-length  cross-section  when  lateral 
gas  flow  is  considered  (i.e.  for  Case  2).  It  is  seen  that  there 
exists  a  large  pressure  drop  in  the  GDL  between  the  adjacent 


Fig.  9.  Pressure  distribution  in  the  cathode  GDL  and  average  velocity  scaled  by  the  cathode  inlet  velocity  at  the  mid-length  cross-section  with  convection  in 
the  GDL. 
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channels  in  different  flow  paths  at  the  two  ends  of  the  fuel  cell. 
Between  Channels  1  and  2,  c.a.  1 .5  kPa  pressure  difference  is 
present.  According  to  the  Darcy’s  law,  Eq.  (5),  and  using  the 
GDL  permeability  of  10“ 12  m2,  one  obtains  a  lateral  velocity 
of  about  0.15  m/s  in  the  GDL  from  the  dry  Channel  1  to  wet 
Channel  2.  As  a  result,  there  is  a  certain  amount  of  inlet 
fresh  air  bypassing  through  the  GDL  and  directly  flowing 
into  the  channel  in  the  other  flow  path  towards  the  exit.  This 
gas  bypass,  or  “short  circuit”,  greatly  reduces  the  reactant 
utilization  in  the  cell. 

To  quantify  the  bypass,  average  gas  velocity  in  the  cathode 
gas  channels  is  also  presented  in  Lig.  9.  It  can  be  seen  that 
bypass  causes  severe  leakage  of  injected  flow.  Lor  Path  1, 
from  the  inlet  Channel  1  to  the  middle  of  the  fuel  cell,  e.g. 
Channel  17,  the  mass  of  flow  in  the  channel  decreases  by 
about  80%.  Prom  the  middle  to  the  outlet  Channel  36,  the 
flow  acquires  four  times  mass  from  the  other  flow  path. 

Besides  severe  gas  bypass,  the  lateral  convection  from  the 
dry  inlet  channel  to  the  wet  outlet  channel  opposes  to  wa¬ 
ter  diffusion,  thus  greatly  diminishing  the  self-humidification 
benefit.  The  effect  of  convection  on  species  transport  is  in¬ 
dicated  in  Pig.  10,  which  superposes  H2O/O2  concentration 
contours  with  the  velocity  vector  plot  in  the  cathode  side  be¬ 


tween  Channels  1  and  2.  It  is  seen  that  about  0.1  m/s  veloc¬ 
ity  in  the  porous  GDL  is  induced  by  the  pressure  difference 
between  the  two  channels,  in  consistency  with  the  estimate 
made  earlier.  If  we  define  the  Peclet  number  as  a  parameter 
to  measure  the  relative  strength  of  convection  to  molecular 
diffusion,  then: 


Pe  =  u 


^land 


D 


eff 

GDL 


Lor  Deff  ~  10-5  m2/s  and  8\md  =  0.001  m,  Pe  is  about  10, 
demonstrating  that  convection  is  dominant  in  the  lateral 
species  transport  through  the  GDL.  As  a  result,  the  bypass 
flow  significantly  offsets  the  internal  humidification  induced 
by  diffusion.  In  addition,  bypass  loses  fresh  air  of  higher  con¬ 
centration  in  O2  to  the  outlet  channel  towards  the  exit  from 
the  fuel  cell. 

Similarly,  strong  bypass  exists  between  other  inlet  and 
outlet  channels,  as  shown  in  Pig.  11  between  Channels  5 
and  6.  The  magnitude  of  the  lateral  velocity  is  reduced  to 
about  0.05  m/s  but  convection  still  dominates  the  species 
transport. 
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Fig.  10.  (a)  H2O  and  (b)  O2  distributions  and  velocity  in  the  cathode  side  at  the  mid-length  cross-section  between  Channels  1  and  2  with  convection  in  the 
GDL  at  Vceii  =  0-65  V  and  /avg  =  0.88  A/cm2. 
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Fig.  11.  (a)  H2O  and  (b)  O2  distributions  and  velocity  in  the  cathode  side  at  the  mid-length  cross-section  between  Channels  5  and  6  with  convection  in  the 
GDL  at  Vceii  =  0-65  V  and  /avg  =  0.88  A/cm2. 


Fig.  12.  Average  H2O/O2  concentrations  in  the  cathode  gas  channel  for  each  flow  path  with  convection  in  the  GDL  at  Vcen  =  0.65  V  and  7av g  =  0.88  A/cm2. 
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Fig.  13.  Water  content  profile  at  the  interface  of  the  membrane  and  cathode  catalyst  layer  at  the  mid-length  cross-section. 


An  integral  effect  of  gas  bypass  is  shown  in  Fig.  12,  which 
plots  H2O/O2  concentrations  in  the  cathode  gas  channel  for 
each  flow  path.  Contrary  to  Case  1 ,  along  each  flow  path,  H2O 
concentration  increases  first  due  to  ORR  water  production 
but  decreases  later  due  to  addition  of  drier  gas  from  the  other 
flow  path.  Furthermore,  water  concentration  reaches  as  high 
as  30  mol/m3  in  the  middle  of  the  fuel  cell,  due  partly  to  the 
deceleration  of  channel  flow  as  shown  in  Fig.  9. 

The  trends  of  O2  concentrations  of  the  two  flow  paths  are 
also  different  from  Fig.  3.  It  is  shown  that  O2  concentration 
decreases  first  along  the  flow  due  to  ORR  consumption,  then 
increases  due  to  bypass  of  fresher  air.  In  addition,  high  O2 
concentration  at  the  outlet  indicates  that  bypass  causes  a  very 
serious  leakage  of  O2  and  hence,  lower  cell  performance  in 
this  case.  Furthermore,  because  O2  concentration  becomes 


very  low,  around  2  mol/m3 ,  in  the  middle  of  the  flow  field, 
O2  transport  limitation  occurs  and  leads  to  low  current  density 
locally. 

The  benefit  of  self-humidification  is  best  reflected  by  the 
membrane  water  content  profile.  Fig.  13  shows  such  profiles 
at  the  interface  between  the  membrane  and  cathode  catalyst 
layer  for  both  cases.  On  the  two  ends  of  the  fuel  cell,  it  can 
be  seen  that  the  membrane  hydration  in  Case  1  is  better  than 
Case  2  where  internal  humidification  is  suppressed  by  flow 
bypass. 

Fig.  14  shows  the  current  density  distributions  at  the  same 
cross-section.  Focusing  on  the  inlet  Channels  1  and  35,  it  can 
be  seen  internal  humidification  increases  the  current  density 
for  Case  1,  comparing  to  Case  2.  In  addition,  bypass  forces 
the  dry  air  across  the  GDL  between  the  adjacent  inlet  and 
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Fig.  14.  Current  density  distribution  profile  at  the  mid-length  cross-section. 
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Fig.  15.  Current  density  contours  for  Case  1  at  Vceii  =  0.65  V  and 
/aVg  =0.91  A/cm2. 

outlet  channels  so  that  lower  current  density,  resulting  from 
higher  ionic  resistance,  appears  under  the  land  in  Case  2. 
As  another  result  of  bypass,  higher  current  densities  under 
Channels  2  and  36  are  seen  in  Case  2  because  02-rich  air  is 
added  into  these  moisture  channels  by  bypass.  Meanwhile, 
there  exist  very  low  current  densities  under  lands  for  Case 
2,  such  as  the  land  between  Channels  14  and  15.  This  is 
because  O2  concentration  reaches  a  low  level,  as  shown  in 
Fig.  12.  Highly  non-uniform  current  distribution  is  apparent 
from  channel  to  channel  in  this  flowfield. 

Fig.  15  shows  the  current  density  contours  for  Case  1. 
This  local  performance  of  the  PEFC  is  mainly  determined 
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Fig.  16.  Current  density  contours  for  Case  2  at  Vceii  =  0.65  V  and 
4vg  =  0.88  A/cm2. 


by  two  factors,  O2  concentration  and  membrane  hydration. 
The  asymmetry  of  the  distribution  between  the  two  counter 
flow  paths  on  the  cathode  is  remarkable,  obviously  due  to  the 
co-flow  configuration  of  the  anode  where  the  anode  inlet  is 
always  located  on  one  side  of  the  anode  plate.  In  addition,  the 
highest  current  density  occurs  near  Channels  2  and  3  because 
these  locations  are  near  the  wet  cathode  outlet  and  anode 
inlet.  Under  the  inlet  channels  in  cathode,  such  as  Channels 
1  and  35,  the  current  density  is  low  due  to  dry  air  injection. 
Along  the  flow  path,  the  current  density  increases  because 
water  production  and  self-humidification  raise  the  membrane 
hydration  level.  This  trend  is  clearly  shown  in  Path  2.  For  Path 
1 ,  more  O2  consumption  near  the  inlet  results  in  decreasing 
current  density  in  the  latter  part  of  this  flow  path.  Fig.  16 
shows  the  current  density  contours  for  Case  2. 

4.  Conclusions 

Numerical  simulations  of  a  double-path  flowfield  PEFC 
with  36  channels  have  been  carried  out  to  study  the  complex 
flow  phenomena  and  water  transport  in  PEFCs.  The  simula¬ 
tion  results  reveal  that  the  effect  of  internal  humidification 
induced  by  diffusion  promoted  by  making  the  two  neighbor¬ 
ing  channels  in  counter  flow  is  significant,  but  is  in  reality 
offset  by  the  opposed  gas  bypass  flow  induced  also  by  the 
same  flow  configuration.  In  addition,  the  bypass  phenomenon 
leads  to  significant  leakage  of  gas,  substantially  reducing  the 
reactant  utilization.  To  decrease  the  effect  of  bypass,  small 
permeability  of  GDLs  was  recommended  for  the  double-path, 
counter-flow  PEFCs. 

While  the  counter  flow  between  two  neighboring  flow 
paths  in  the  same  cathode  flowfield  is  found  ineffective 
through  this  extensive  numerical  study,  this  configuration 
provides  a  unique  opportunity  to  explore  the  rich  physics 
behind  the  complex  flow  and  transport  phenomena  occurring 
in  a  PEFC  as  well  as  to  test  any  numerical  models  of  PEFCs. 
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